Direct Numerical Confirmation of Pinning Induced Sign Change in Superconducting 

Hall Effect in Type-II Superconductors 
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Using the time-dependent Ginzburg-Landau equation with the complex relaxation time and the 
Maxwell equation, we systematically examine transverse motion of vortex dynamics in the presence 
of pinning disorders. Consequently, in a plastic flow phase in which moving and pinned vortices 
coexist, we find that the Hall voltage generally changes its sign. The origin of the sign change is 
ascribed to a fact that moving vortices are strongly drifted by circular current of pinned vortices 
and the enforced transverse moving direction becomes opposite to that by transport current. This 
suggests that the Hall sign change is a behavior common in all disordered type-II superconductors. 
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Since the discovery of cuprate high-T c superconduc- 
tors, much attention has been devoted to vortex dynam- 
ics in not only superconductors but also various superflu- 
ids from liquid Helium to atomic gas. In particular, vor- 
tex pinning dynamics under disorders inevitable in super- 
conducting materials is a central issue of vortex physics 
because of deep relations to its industrial applications. 
In this paper, we numerically examine the vortex pinning 
dynamics and give a clear explanation to a controversial 
topic in vortex physics, i.e., sign change in the supercon- 
ducting Hall effect 0, S, S, 0, H . 

The equation of motion for a moving vortex in super- 
fluids has been highly controversial [5j . The heart of the 
problem is the non-dissipative transverse force (or the 
vortex velocity part of Magnus force) 0,0,110, which 
brings about the superconducting Hall effect. In BCS su- 
pcrfluids, the transverse component is generally subtle, 
and therefore, various theoretical proposals remain un- 
examined fully. Among those theories, one of the most 
controversial struggles is whether pinning or disorder can 
be an origin of the sign change in the superconducting 



Hall effect [13 El El, E3, El EE El] . If it is true, it then 



indicates that the sign reversal is not limited in particular 
superconductors, but universal for all disordered type-II 
superconductors. In this paper, we clarify that the idea 
is really true by numerically solving the time-dependent 
Ginzburg-Landau (TDGL) equation with complex relax- 
ation time and the Maxwell equation. This is a direct 
and clear confirmation of the pinning induced sign change 
without simplification and modeling. 

The vortex dynamical phases under disorders are 
roughly classified into two types, i.e., plastic and col- 
lective flow phases [l7|. In the former phase, moving 
and pinned stationary vortices coexist. Then, it is found 
that moving vortices are strongly drifted by the circu- 
lar current of pinned vortices and its inducing transverse 
moving-direction can become opposite to that by trans- 



port current. Such a feature is schematically displayed in 
Fig. 1(a). Thus, the averaged transverse moving direc- 
tion in the plastic flow phase can be different from that 
in the collective one. This is an origin of the Hall sign 
change confirmed in this paper. 

Let us present the system setup to confirm the pinning 
induced sign change. We prepare a two dimensional sys- 
tem in the xy plane. The external magnetic field H a is 
applied perpendicular to the plane. To simulate vortex 
dynamics, we numerically solve the the TDGL equation 
coupled with the Maxwell one written as [18|, 
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Here, we introduce local suppressions of the transition 
temperature T c (r), which act as vortex pinning sites. 
The order parameter A is normalized by its mean field 
value at the zero temperature without the magnetic field, 
Ao, and time t, vector potential A, and magnetic field 
H are done by to — <in n 2 / c 2 , A a — o /(27r£ o ), 
and Ho = 0o/(27r£ 2 ), respectively, where £o, K , °", c, 
and <f>o{= lithe/ (2e)) are the zero-temperature coherence 
length, the Ginzburg-Landau parameter, the normal- 
state longitudinal conductance, the light velocity, and 
the flux quantum, respectively. To keep the gauge invari- 
ance of Eqs. (1) and (2) on numerical grids, we use the 

link variable U l J = exp 

stands for x or y 1 
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. The magnetic field H 
is given by the Stokes' theorem J S (H/H ) ■ ndS/H = 
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FIG. 1: (Color online) (a) A schematic figure of a moving single vortex driven by an applied current in the case of a complex 
TDGL relaxation rate in the presence of the pinned vortex, (b) The contour map of the local transition temperature T c (r). 
The closed curves (black) show the positions of pinning sites, and the shading color indicates degree of the T c suppression, (c) 
The superimposed snapshots of moving vortices in the time interval 4 x 10 s < t/to < 6 x 10 s . The vortices are visualized at 
each instant by the contour plot (blue) of the order parameter at A/Ao = 0.2. The positions of pinning sites are marked by 
black curves. T/T c0 = 0.69, H a /H = 0.2, and j x /j = 5 x 10~ 5 . 




FIG. 2: Temperature dependences of the longitudinal (V x ) and Hall (V y ) voltages in units of Ao^o/to- The filled (open) circles 
indicate the longitudinal (Hall) voltage. Note that the Hall voltage is plotted in the ten times larger scale. The voltages are 
calculated from the electric field averaged over —85 < r x {y)/£,a < 85 and 4 x 10° < t/to < 6 x 10°. The applied current is 
jx/jo = 5 x 10 -5 . The applied field is H a /H = 0.15 for (a) and (c), and H a /H = 0.2 for (b) and (d). The results in (a) and 
(b) are obtained for the distribution of pinning sites shown in Fig. 1(b), while the same distribution is used but whose bottom 
is inverted into top for moving vortices to obtain the results (c) and (d) (see text). 



J (A/Aq) ■ dl/^o, and the electric field is calculated by 
i P = ~(Ao/t )f s d(A f t/A )/d(t/t Q )d 2 r/S, where S is 
the unit plaquette surrounded by link variables. We eval- 
uate the longitudinal and the Hall voltage from E^. In 
order to concentrate on the vortex contribution to the 
Hall voltage, we neglect the normal-state Hall conduc- 
tivity in Eq. @ for clarity. Instead, the dimensionless 
relaxation rate Y in Eq. |T]) is set to a complex number. 
According to Ref. [2^ |. \/Y is related to the forces act- 
ing on each moving vortex. If we set 1/T pure real, the 
transverse force and the resulting Hall voltage are zero. 
On the other hand, a finite imaginary part of 1/T brings 
about a transverse force, and the sign of the im aginary 
part controls the sign of the Hall effect @, HI 0, [2J, Ej- 
We keep the imaginary part positive and never change its 
value throughout this study, i.e., the transverse force is al- 



ways positive and unchanged. Under the condition of the 
fixed complex relaxation rate, we find that the Hall volt- 
age amplitude diminishes and a sign reversal occurs in 
the plastic flow regime owing to moving vortices affected 
by the circular current around pinned vortices. 

In the present simulation, the system size is 200£o x 
200£o, which is discretized by the square grid whose unit 
dimension is £o x £o ■ The external current is applied along 
the x direction, and a periodic boundary condition is im- 
posed in this direction to eliminate edge boundary effects 
on the vortex motion. The remaining boundary edge per- 
pendicular to the y direction relevant to vortex entry and 
escape is modeled as an interface between a superconduc- 
tor and a normal metal. Around this interface, T c (r) is 
set T c /T c q = 0.1r/£o) where r is the distance from the 
interface and T c q is the bulk value of T c . The number 
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FIG. 3: (Color online) (a) The distribution of pinning sites in the focused area, whose location is marked by the small box 
in Fig. 1(b). (b)The transverse electric field in the area for T/T c = 0.69, H a /H = 0.2, and j x /jo = 5 x 10" 5 . The data is 
averaged over 5.8 x 10 5 < t/to < 5.84 x 10 s . (c)The vortex positions at t/to = 5.8 x 10 5 and 5.84 x 10 5 through the contour 
plot of the order parameter amplitude. The arrows indicate the moving directions (d)The schematic representation of moving 
vortices around stationary vortices. 



of vortex pinning centers is 500 inside the present system 
200£o x 200£o- The size of each pinning center is 2£o x 2£o, 
inside which T c (r) is randomly suppressed in the range 
0.8 < T c /T c q < 1. The locations of pinning centers are 
randomly distributed, e.g., as shown in Fig. 1(b). 

In the TDGL dynamical simulation, we prepare an ini- 
tial state in the absence of both the applied current and 
external field, and then start to apply a current j x and 
a target external field H a at t = 0. The applied cur- 
rent density is always set as jx/jo = 5 x 10 _J , where 
jo(= 0o/(27t£q)) is the depairing current. The GL param- 
eter is k — 2.83, and the minimal time step is 3 x 10 _3 io- 
We fix 1/r = 1 + 0.3i, which leads to a substantial ratio 
of the Hall voltage V y to the longitudinal one V x , i.e., 
Vy/V x ~ 0.2, in the uniform current under no pinnings. 
Such a large imaginary value can give a striking contrast 
in the sign change of the Hall voltage. To avoid counting 
an interface influence on the voltage, e.g., an effect of the 
diamagnetic current, we take an average of the electric 
field within the region —85 < r x r y )/£o < 85. The time 
average of the Hall voltage is taken over the time interval 
4 x 10 5 < t/to < 6 x 10 5 , during which vortex motions 
are steady. 

Let us present simulation results. First, the tempera- 
ture (T) dependences of the longitudinal and Hall volt- 
ages under the applied fields H a /H = 0.15 and 0.2 are 
displayed in Figs. 2(a) and 2(b), respectively. The longi- 
tudinal voltage V x monotonically decreases with decreas- 
ing T in both cases. Although the Hall voltage V y ex- 
hibits non-monotonic behavior, both the signs become 
negative in the region of the small longitudinal voltage 
V x . Here, one might imagine that this negative transverse 
voltage occurs because of guided vortex flow lines [l3| 
formed accidentally by clustering of pinned sites. How- 



TABLE I: Yes/No table in terms of the sign reversal of the 
Hall voltage for different random pinning distributions (I)- 
(IV) and their field inversion cases (I)— (IV). "Y(N)" means 
that the sign reversal is observable (or not). 
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ever, it is not always the case. In order to confirm it, we 
repeat the simulation by just reversing the magnetic field 
direction only. While the sign reversal disappears for the 
field H a /Ho = 0.15 (Fig. 2(c)), it is kept for another 
field HJHq = 0.2 (Fig. 2(d)). This suggests that there 
must be an intrinsic origin of the sign reversal beyond 
the guiding effects. We further perform simulations (not 
shown) under other random vortex-pinning distributions 
(II)-(IV) and repeat the simulations with the field direc- 
tion reversal (II)-(IV). These results are summarized in 
Table I. The guiding vortex flow state is found to indeed 
cause the sign change under the distributions II and III, 
because the Hall sign does not change on the field inver- 
sion cases in both II and III. However, the sign reversal 
is kept on the field direction reversal for the other dis- 
tributions I and IV. From these results, it is deduced 
that the intrinsic sign reversal effect always coexists with 
the guiding one in the plastic flow phase. One naturally 
expects that the intrinsic mechanism dominates in suffi- 
ciently large samples. The following analysis attributes 
the intrinsic mechanism to an effect of the circular cur- 
rents around pinned vortices. 

Let us focus on detailed vortex dynamics around 



4 



vortex-pinning sites. Figure 3(a) is an enlarged figure 
of a small area whose location is marked in Fig. 1(b). 
The distribution of the transverse electric field averaged 
over 5.8 x 10 5 < t/t Q < 5.84 x 10 5 is shown in Fig. 3(b), 
where the sign of the transverse field is negative near 
the vortex pining sites. In addition, by monitoring the 
vortex position at 5.8 x 10 5 to and 5.84 x 10 5 io in the dis- 
played area, the focused vortices are found to flow against 
the applied transport current [Fig. 3(c)]. These observed 
results signify the following. The pinned vortex has a 
circular current flow around its core. The flow direction 
is opposite to that of the transport current in its top half 
as schematically shown in Fig. 3(d), where locally the 
Lorentz force (or the superfluid velocity part of Magnus 
force) is directed from bottom to top. Then, the moving 
vortex penetrating into the circular-current flow range is 
drifted into the opposite direction to the transport cur- 
rent. This is because the positive imaginary part of T 
fixed in this paper always drives the moving vortex into 
the downstream side of the current flow as usual. Thus, 
the moving vortex exhibits a motion directed opposite to 
the transport current. 

Finally, let us discuss the present results through a 
comparison with experiments. The present calculations 
have revealed that when the vortex dynamics change 
from the collective flow to the plastic flow phase the mov- 
ing vortex frequently reverses its transverse moving di- 
rection. This directional reversal principally requires the 
plastic flow phase as a vortex dynamical phase. In other 
words, this effect is universal for all type-II superconduc- 
tors as long as disorders or pinning sites enough to keep 
the plastic flow phase are introduced inside the sample. 
In high-T c cuprate superconductors, double sign change 
in addition to single one have been frequently observed 
depending on the sample 26] . These experimental results 
can be explained on the basis of the present result as fol- 
lows. When the current carrying phase changes from the 
normal to the flux flow phase, the first sign change oc- 
curs. This can be interpreted by an idea that there is 
a difference between the Hall effect in the normal phase 
and the fluctuation Hall effect near the superconduct- 
ing transition via the relaxation of the order parameter 
@j 0, 23, 24, (2||. This is a microscopic sign change 
mechanism depending on the electronic structure. On 
the other hand, it has been observed that the final rever- 
sals strongly depend on the sample quality or the rate 
of artificial damage to enrich pinning centers [271 ]. Thus, 
the final sign change is attributed to the pinning induced 
one as confirmed by the present simulation. Moreover, 
such a sign change is also well-known to be sensitively 
dependent on the sample quality in conventional type-II 
superconductors. 



In conclusion, we performed TDGL simulations with 
the complex relaxation time to confirm the pinning 
induced sign reversal of the superconducting Hall ef- 
fect. Consequently, the simulation revealed that the sign 
change can occur when the current carrying state enters 
the plastic flow phase from the collective flux flow one. 
Moreover, the detailed analysis on the vortex motions 
successfully explained that, when the circular current of 
the pinned vortex strongly drifts the moving vortex in 
the plastic flow phase, the moving vortex feels the trans- 
verse force causing the Hall sign change. These results 
suggest that the Hall sign change is an indicator of vortex 
dynamical phases in disordered type-II superconductors. 
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See Table I in Ref. [IJ . Two limiting cases of the dirty and 
clean does not show the sign change of the Hall voltage 
in V and Nb. 



